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Abstract. Minimal surfaces in the sub-Rieniannian Heisenberg group can be constructed by means of 
a Riemannian approximation scheme, as hmit of Riemannian minimal surfaces. We study the regularity 
of Lipschitz, non-characteristic minimal surfaces which arise as such limits. Our main results are a-priori 
estimates on the solutions of the approximating Riemannian PDE and the ensuing C°° regularity of the 
sub- Riemannian minimal surface along its Legendrian foliation. 



1. Introduction 



The first Heisenberg group H^ is a Lie group with a 3-dimensional Lie algebra h — V^ (B V^ such that 

dim(yi) = 2,dim(y2) = i, [V^^V^] = V^ and [h,V^] = 0. Let S,X,y £ h he any basis such that 

[5, A"] = y E V^. By assigning a left-invariant Riemannian metric go on the horizontal sub-bundle HM.^ 

given by the V^ layer, we obtain a sub-Riemannian space {M^,go). We choose S,X such that they are 

orthonormal with respect to gg. The corresponding control metric dg (the Carnot-Caratheodory metric j32l) 

K^ I is easily shown to be well defined. We extend go to a (left-invariant) Riemannian metric gi on the full tangent 

^y ■ bundle of h requiring that V^ and V^ are orthogonal in this extension. The dilated metrics g^, £ > are 

defined so that S, X, ey are orthonormal. We define d^ to be the corresponding distance function. We define 

^■f-N ■ polarized coordinates {xi,X2,X3) in H^ by identifying the triplet with the point e'Kp{x3S) eicp{xiX + X2y)- 

.' I The Baker-Campbell-Hausdorff formula yields S = ds, X = di + X3d2 and y = 82- 

(<— ^ ■ If M C H^ is a C^ surface, then p & M is called characteristic if both 5, X are tangent to M at p. An 

QQ ' intrinsic graph (see [HI US]) is a (non-characteristic) graph of the form 



(LI) M - {.T3 - u{xi,X2)\ (xi,X2) enc M^}, 



^^ . where fi C R^ is an open set. An analogue of the classical implicit function theorem [25l [HI [T] shows that 

^ I any surface {/ = 0} with Sf, Xf e C(IHI^) can be represented as an intrisic graph, in a neighborhood of any 

of its non-characteristic points. 

The flow associated to the line bundle of tangent directions which are also horizontal foliate the com- 
plement of the characteristic locus of the surface. This is called Legendrian foliation in the literature. 
Note that the horizontal tangent bundle of an intrisic graph (jl.ip is spanned by the single vector field 
T = XuS -\- '^|„((xi,x2),z:i.a;2)- ^c note that the projection of this vector field on Til yields the vector field 
in f7, 

Xi^u =di +ud2. 
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Minimal surfaces. Several equivalent notions of horizontal mean curvature Hq for a C^ surface M C H^ 
(outside characteristic points) have been given in the literature. To quote a few: Hq can be defined in terms 
of the first variation of the area functional [5D1 [221 IHl [3S1 113 [SB IZ] , as horizontal divergence of the horizontal 
unit normal or as limit of the mean curvatures H^ in the Riemannian metrics g^ as e ^> 0. It is also well 
known (see for example [9l[29l[T8]) that i/o coincides with the curvature of the projection of the Legendrian 
leaves on the Horizontal plane. 

A C^ non characteristic surface M C H^ is called minimal if it satisfies Hq = identically. In particular 
for a C^ intrinsic graph, a direct computation yields that the PDE can be written in terms of the vector 
Xi_u as follows 

(1.2) gp ^ Xi ^ I ^ ^^-"" 1^0. 

A deep result of Ambrosio, Serra-Cassano and Vittone IJ shows that such PDE continues to hold below the 
C^ threshold in a suitably weak sense. 

Generalized solutions and Riemannian approximants. Because the horizontal mean curvature arises 
as first variation of the sub-Riemannian perimeter, minimal surfaces are critical points of the perimeter. As 
such these objects can be interpreted in weak sense, far below the threshold of C^ smoothness (see [T], [35], 
[231j[211j[9],[II],[35]). As an example, starting from the family of shears X2 — xix^ +5(2^3) (see [34]) one can 
obtain the intrinsic graph ^3 = u{xi,X2) = ^ -sqn(x ) defined in Jl = {(a;i,a;2) G M^ : xi > 1} which is 
minimal in the sense that it is foliated by horizontal lifts of segments, it is locally Lipschitz (with respect to 
the Euclidean metric) but clearly not C^ smooth. 

In [33] and in [TT], the authors prove existence of (respectively W^'^ and Lipschitz) minimal surfaces using 
the Riemannian approximation scheme: as e -^ 

(H\4)^(H\do), 

in the Gromov-Hausdorff topology (see [71 Section 2.4] for a detailed description). 

The approach to existence of solutions in these papers is based on a-priori estimates for the minimizers 
of the approximating Riemannian g^ perimeter functionals [33] and on solutions of the ^''Riemannian" regu- 
larized versions of (|1.2p [11] . In adapting the approximation scheme to the intrisic graphs setting we note 
that the minimal surface PDE for the metric g^ corresponding to intrinsic graphs (|l.ip 



L.u.±X..^- 



(1.3) L,u = ^x,„ ( ,/';"! \ = 0, in n C 



where 

A"2,u = £^2 Ve = {Xi,u,X2^u), 

is a natural elliptic regularization of the PDE (|1.2[) . As such, it is more amenable to establishing a-priori 
higher regularity estimates. The difficulty of course resides in obtaining estimates which are uniform in the 
parameter e as e ^ 0. 

These observations lead us to the definition of the class of minimal surfaces we want to investigate 



Definition 1.1. We say that a function u S LipiVL) is a vanishing viscosity solution of the equation M.2fl 
if there exists a sequence of positive numbers Ej with Sj — > when j —^ c», and a sequence (uj) in C°°{Q) 
such that: 

(V Ell X^.u, ( ^^v:"l,p ) =0. ^nn for all j G N. 

(a) The sequence (uj) is bounded in Lip{Q) and uniformly convergent on subcompacts of fl to u. 
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Remark 1.2. Existence of vanishing viscosity solutions in the case of t— graphs, i.e. graphs of the form 
X2 — g{s,xi), has been proved in |1H Theorem A and Theorem 4.5]. In the same paper the authors establish 
that such solutions are perimeter minimizers and address uniqueness. Assuming C^ convergence of the 
approximating solutions, outside the characteristic sets the t— graphs solutions in jllj , both the approximating 
and the limit solutions, can be represented as intrisic graphs and hence yield vanishing viscosity intrinsic 
graphs. 

Regularity results. Given the examples of non-smooth minimal surfaces mentioned above, the question 
arises as which kind of regularity can one expect. This problem has beeen recently addressed in a series of 
papers [34] , [11] , [5] , [TO] and [3] . Valuable insights into the problem of regularity are also provided in the 
works j34j . [llj . and |35j in the form of examples of non-smooth minimal surfaces. The regularity properties 
of the implicit function in the implicit function theorem quoted earlier provides an interesting insight into this 
problem and indicates that one should look for regularity only in the direction of the Legendrian foliation. 
Indeed we prove 

Theorem 1.3 (Regularity). If u is a vanishing viscosity solution of (|1.2[) . then for all a G (0,1) and 

K ccnc R^, 

(1.4) ueC^'^iK) 

and for all k G N and p > 1 

(1-5) xl^uewl^^m. 

Here C^'"{K) and Wi^^{fl) denote the spaces of functions with Holder continuous Euclidean gradient and the 
classical Sobolev Space. In particular one has X^ ^u G Cj^^{n) (the Euclidean Holder space) for all a G (0, 1) 
and hence X^ ^u = 0, holds pointwise everywhere. 

Remark 1.4. To better understand the notion of intrinsic regularity we return to the non-smooth minimal 
graph u{xi,X2) = -s^n (x ) described earlier. Although this function is not C^ in the Euclidean sense, 
observe that Xiu = for every a:i,a;2 G fi Hence, this is an example of a minimal surface which is not 
smooth but which can be differentiated indefinitely in the direction of the Legendrian foliation. 

Remark 1.5. The regularity theory for intrinsic minimal surfaces in H" with n > 1 is quite different. In 
the recent paper 5 we show that any Lipschitz continuous vanishing viscosity minimal intrinsic graph in 
H", n > 1 (defined through the Riemannian approximation scheme) is smooth. The main reason is that in 
higher dimension the horizontal tangent bundle generates as a Lie algebra the full tangent bundle, while this 
does not happen in the n = 1 case. 

As a consequence of the regularity theorem we can prove that the Sobolev weak derivatives of vanishing 
viscosity solutions agree with Lie derivatives along the leaves of the Legendrian foliation. Hence, we obtain 
that vanishing viscosity solutions actually satisfy (|1.2p everywhere pointwise. This result immediately yields 
a rigidity of the Legendrian foliation. 

Corollary 1.6 (Lie differentiability and Legendrian foliation). Letx^ — u{xi,X2), (xi,X2) G fJ 6e a Lipschitz 
continuous vanishing viscosity minimal graph. The flow of the vector Xiu yields a foliation of the domain 
f2 by polynomial curves 7 of degree two. For every fixed xq G fi denote by 7 the unique leaf passing through 
that fixed point. The function u is differentiable at xq in the Lie sense along 7 and the equation ()1.2|) reduces 
to 
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Comparison with other regularity results. We describe the relation between our results in Theorem 
11.31 and Corollary 11.61 and the regularity results in [10] and in [3] . In [10] Cheng, Hwang and Yang prove 
that any C^ weak solution of the prescribed (continuous) horizontal mean curvature PDE, has C^ smooth 
Legendrian foliation outside of the characteristic set. In ^, Bigolin and Serra Cassano study the regularity 
of minimal intrisic graphs (jl.l[) where (jl.2[) is interpreted in a weak sense (i.e. broad* solutions defined in 
|3l Definition 3.1]) and prove (among other results) that Lipschitz regularity of the intrisic gradient Xi^^u 
implies the Euclidean Lipschitz regularity of the function u. In the present paper we require only Lipschitz 
continuity of m and prove higher order intrinsic differentiability than either [3], and [TU]- On the other hand, 
we only deal with the case Hq = and with those solutions which are limits of Riemannian minimal graphs. 
In this sense our results are more specialized than the ones in the other two papers. 

Applications of Theorem ll.3l Invoking the implicit function theorem, we can apply Theorem ll.3l to study 
the regularity away from the characteristic locus of the Lipschitz perimeter minimizers found in [111 . Since 
the results in that paper apply to t — graph and not intrinsic graphs we need some extra assumptions on 
the convergence of the approximating solution to be able to invoke our intrinsic graphs regularity. Here and 
in the following Vb denotes the Euclidean gradient in R^, and (21,2:2,23) are the exponential coordinates 
exp(znS + Z2X + z^y) = exp(2;35) exp(xiA' + X23^). 

Corollary 1.7. Let O CZ M? be a strictly convex, smooth open set, (f> e C^'"(0) and for each (21,2:2) G O 
denote by (21, 22)* — (22,-21). Consider the family 

{9s{zi, Z2)}e sup \ge\ + sup \V Ege\ < C (uniformly in e), 
o o 

of smooth solutions of the approximating minimal surface PDE 

,. ( Vis.ge + (21,22)* \ n ■ r, .1 A, ■ -^n 

div — 7^^^^^^^^^^^^^^= = V m (J and g,, — (p in o(J 

\^e^ + \yEge + izi,z2)*\J 

found in 11, Theorem 4.5]. If for po = {pI,Pq) E O, a > and for every e > we have \dzige{po)\ > a > 
(or any other partial derivative is non-vanishing at po uniformly in e) then there is a sequence e^ — > 
such that the Lipschitz perimeter minimizer g = liuiek^o 9sk satisfies g € C^'°' and is infinitely many times 
differentiable in the direction of the Legendrian foliation of z^ — g{zi, 22), "in a neighborhood of the point po- 

Proof. The implicit function theorem implies that the level set of 

23 -56(^1,^2) 

can be written as smooth intrinsic graphs a;3 — Ue{xi,X2) in a neighborhood fJ of (pj,pg, (7(730) )• The Lipschitz 
bounds on g^ (proved in 11, Propositions 4.2-4]) yield uniform Lipschitz bounds on u^, thus allowing to 
apply Theorem 11.31 and conclude the proof. D 

From this result one may conclude immediately the following 

Corollary 1.8. Let 23 = (7(21,22), (21,22) € f2 C M^ be a C^ minimal graph which is the C^ limit of 
Riemannian minimal graphs as in \\\\ Theorem 4.5]. In the neighborhood of any non- characteristic point 
g € C^'" and is infinitely many times differentiable along the Legendrian foliation. 

Theorem 11.31 can be also used to rule out minimal intrinsic graphs which do not arise as limits of Rie- 
mannian minimal graphs. For instance, the example exhibited above 0:3 = u{xi,X2) — — _ ^^ / — y defined in 
Vt = {(xi, 0:2) S M^ : xi > 1} cannot be a vanishing viscosity minimal graph as it lacks the C^'" regularity 
from Theorem ll.3l 
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Sketch of the proof and final remarks. We now turn to a description of the techniques used in the 
proof of Theorem 11.31 Since u is a vanishing viscosity solution, this theorem is proved by means of a priori 
estimates (uniform in the parameter e as it decreases to zero) for each element of the approximating sequence 
Uj of solutions of (II. 3p . 

The PDE (|1.2p . has a structure similar to the Levi equation in M.^ . In fact, the Levi equation can be 
represented as Riemannian approximation of a sum of squares of vector fields 

Zl + Wuu^ + e^dlu ^ 0, 

for suitable non linear vector fields Z„, Wu depending on the solution u. Regularity results for Lipschitz 
continuous viscosity solutions were established in [13 and [16j. The techniques used in these papers are 
based on a modification of the Moser iteration, along with representation formula and uniform estimates 
on the fundamental solution. In [16j . [19j the authors address properties of the analogue of the Legendrian 
foliation for a Levi flat graph. 

The cited work on the Levi equations provides a coarse outline and a strategy for the proof of Theorem 
11.31 However, the equation ()1.2p presents additional difficulties (lack of a background sub-Laplacian, worse 
nonlinearity) with respect to the Levi equation, and the adaptation of the known techniques is very non 
trivial and requires new ideas. 

The first step in the proof involves the study of a linearization (of sorts) of (|1.3p (for simplicity we will 
refer to the approximating functions Uj simply as u) 



(1-6) M,^,,z = Y, X,..u ,^ ^:" ,, = 0, in f7 C 



There are two main difficulties in establishing a-priori estimates, uniform in e, for this PDE: 

• The first problem is due to the fact that the coefficients of the equation involve the function u which 
although smooth satisfies a-priori bounds which are uniform in e only for the Lipschitz norm. To 
deal with this obstacle we operate a freezing argument, substituting the function u with an analogue 
of its ffist order Taylor polynomial, and then carefully study the remainder terms. The regularity 
of rough coefficients degenerate elliptic PDE has been studied by many authors, see for instance the 
monograph [39 for a survey of the literature and new, ground-breaking techniques. 

• The second difficulty stems from the fact that, although ()1.6p is elliptic, its coerciveness degenerates 
as e — > 0. Now, the approximation of degenerate elliptic operators with elliptic regularization 
is a well known and widely used trick. For instance in |30| . the sub-Laplacian Cu — Y^^iLi-^f'^ 
associated to a system of Hormander vector fields is approximated by C^u = Cu + eAu. While it 
is true that the ellipticity of C^ degenerates as e ^ 0, and hence the constants involved in elliptic 
estimates degenerate as well, the operator C/^ satisfies sub-elliptic estimates, which are uniform in s. 
In our case however, the left-hand side of equation (|1.6p approximates not a sub-Laplacian but the 
operator Xi^u{Xi^uz/ y^l + (Xi^„m)^) which is not sub-elliptic. To solve this problem, and obtain 
the regularity in L^ of the derivatives of z, we introduce a completely new ad- hoc lifting process, 
inspired in part by Rothschild and Stein's techniques in [35]. The vector fields Xi_.u,X2,u are lifted 
to a three-dimensional space il x (—1,1) by adding a new variable s and horizontal vector field dg- 
The lifted vectors are Xi = d^i + {u{x) + s'^)dx2, X2 = edx^, and X^ = dg and the set {Xi.X^} 
form a step-three bracket generating system whose commutators yield the direction of degeneracy 
dx2 of (|1.3p . At this point we operate a freezing argument and approximate the operator M^^u with 
higher dimensional Hormander type sub-Laplacians (namely (|3.16p and (|3.17p . built from the frozen, 
lifted vector fields. In this way we extablish a priori W'^'^ estimates uniform in e (Theorem 13. ip . 
which will be the starting point of the regularity proof. 
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In the proof of Theorem 11.31 wc wiU switch back and forth from representations of ()1.3|1 (and its differ- 
entiated versions) in divergence and in non-divergence form. The former works best to deal with higher 
regularity, via Caccioppoh estimates. The second is tailor-made for the freezing technique and the W^'^ 
estimates. 

An important ingredient in the proof is the recent result in [15] . establishing uniform estimates on the 
fundamental solutions of Riemannian regularizations of sub-Laplacians (see the statement in Theorem 13. 3|) . 



For other aspects of minimal surfaces in the Heisenberg group, including classification and Bernstein-type 
results, see [3 [27l El [22l EU [37l [36l [35]. These works also contain more comprehensive lists of references. 

Acknowledgements. The authors would like to thank Manuel Ritore, Francesco Serra Cassano and Paul Yang 
for sharing with them their preprints [35j . [lOj and [3]. 

2. Preliminaries 

In this section we will always assume that u and / are fixed smooth functions defined on an open set fl 
of ]R^,and that m is a solution of the PDE L^m = / in fi. In particular we remark that 

(2.1) IkllL-(n) + l|VeU||Loo(o) + ||S2u||l-(0) < CX), 

and we set 

(2.2) M = ||w||L~(n) + ||Veu||ioo(o) -I- P2u||loo(o) 

where for any function (jj defined on fl we have let Ve(/) = (Xi,„0, X2,u4')- 

We will use the notation W^'P{il.), p > 1 to denote the Sobolev space corresponding to the norm 
||(/)||^^i,PjQs = ||</'||lp(o) + ll'^e0llLp(si)- For simplicity, unless we want to stress the dependence on m, we 

will simply write Xi,X2 instead of Xiu,X2.u- We will denote by Wq'^(17) the space of L'P{n) functions cj) 
such that X[cj) G LP{n) for all 1 < / < fc. 

We recall that, under assumption (|2.ip . the following result holds, (see [S]) 

Proposition 2.1. Let u be a solution of equation il.3]) satisfying 112.1]) . For every compact set K <Z<Z fi 
then there exist a real number a and a constant C , only dependent on the bounds on the constant M in i2. ^|) 
and on the choice of the compact set K such that 

\Mw!-\k) + P2u|lvy^i.2(K) + l|w|lci-°(K) ^ ^■ 
2.1. An interpolation inequality. 

Proposition 2.2. For every p > S, there exists a constant C , dependent on p, and the constant M in 
such that for every function z G C°°(r2) and for every cj) G C^i^D.), and every 5 > Q 

|X,z|P+l/2<^2p ^C I |^|^|4p+2^2p ^ |^|(2p+l)/2|^^^|(2p+l)/2 ^ |^| (2p+l)/2^2p j 
+<S||V,(|X,Z|(^-1)/2)|202P^ 

where i can be either 1 or 2. 
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Proof. This is a slight variant of [131 Prop. 4.2]. We have 

(integrating by parts and using the fact that X^ = —Xi — d2U and X2 — —X2) 

= - f Snd2Uz\X,z\P-'/\ign{X,z)cj)^P - (p- 1/2) f zXfz\X,z\P-^/^<j)^P 



(2.3) -2p f z\X,z\P-^^higTiiX,z)^''P-'X,(f> < 

(where Sij denotes the Kroeneker's delta, by Holder inequality) 

+ Sl |X,zr+l/2^2p + ^1 |V,(|X,Z|(P-1)/2)P02P^ 

choosing 6 > sufficiently small we conclude the proof. D 

A slight modification of the previous proposition, is the following: 

Proposition 2.3. For every p > 3, for every function z € C°°{il,) there exists a constant C , dependent on 
p, the constant M in i2.^) such that and for every (j) S C^ (ft) , and every S > 

\X,z\P+^cb^P < 
<ci I (zP+'6'P + z'\X,zr'cb^P-'\XMA + / \Xfz\'\X,z\P-'\z\'6'P 



where i can be either I or 2. 
Proof. We have 

/ |X,2|P+V'P = I X,z\X,z\Psign{X,z)cf'P = 
(integrating by parts and using the fact that X^ = —Xi — d2U and X| = —X2) 

= - j 5ad2Uz\X,z\Psign{X,z)4>^P~p I zXfz\X,z\P-^(jf'P 

(2.4) -2p I z\X,z\Psign[X,z)4>''P-^X,(j) < 

(by Holder inequality) 

-J J i^"^''^''' + AXz\P-'(b'P-^\X,(b\') +5 j\X,z\P+^(jf'P + jj \z\^\Xfz\'\XMP-^^P, 
choosing 6 sufficiently small we obtain the desired inequality. D 
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2.2. The horizontal mean curvature as a divergence form operator. We now prove that if m is a 
smooth solution of equation p.3|) then its derivatives d2U and Xku are solution of a similar mean curvature 
equation with different right hand side (see also [51 Lemma 3.1]). Differentiating the equation L^u = with 
respect to X^ one obtains 

Lemma 2.4. If u is a smooth solution of L^u — then z — X^u with k < 2 is a solution of the equation 



(2.5) xJ^mB^X.z 



= -{X.M ''''' 



v/l + |Ve«P 






where aij are defined as 

(2.6) fly ■.R'^R a,, (p) = % - Y^^- 

Lemma 2.5. If u is a smooth solution of L^u — then v — d2U is a solution of the equation 



(2.7) XA^^^B^X.V 



aii(Veu) ^,3 aij(Veu) _^^ ^ 



v/l + |VeWP v/l + |VeUP 



aa(V^^, 






Vi + |Vs"l 



2 



where aij are defined in 

Proof. Differentiating the PDE we obtain 

^^[d2.X,](^I^^)-xJ^^lBd=[d„X,]u 

^^,dJ^^^]-x,(^:^Ed=v' 



Vl + |VeUp/ V^l + |V,u| 



D 



Let us consider the linear equation satisfied by the components of the gradient of u: 

M.z^xJ °'^^--^ _ x,z] 



aij{\Ieu) 



If z is a smooth solution of equation 

(2.8) M,z = /, 

then its intrinsic derivatives XiZ are still solutions of the the same equation with a different right-hand side. 
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Lemma 2.6. If z is a smooth solution of 112. 8]) then si — Xiz is a solution of the equation 
(2.9) M,s, = X,f + X,{ °'^(^^") d^uX^)] 

\^/i + \v^uY j 

An analogous computation ensures that 
Lemma 2.7. If z is a solution of \2.8\) then S2 = X2Z is a solution of the equation 



(2.10) M,S2 = X^f X, ( -^B^I^d^uX^z] 

3. The horizontal mean curvature as a non-divergence form operator: Wloc a priori 

ESTIMATES. 

The operator Lg defined in (|1.3p can be represented in non-divergence form 

2 

(3.1) N,u^ Y^ a^j{Veu)X,XjU, 

where Oij are defined in (|2.6p . 

Following the approach in the papers [131 I16j we linearize the operator A^g in the following way: While 
the coefficients of the vector fields Xi depend on a fixed function u, they will be applied to an arbitrary 
function z, sufficiently regular. The associated linear non divergence form operator is 

2 

(3.2) Ni;^uZ^ y^ a.ij{Veu)Xi^uXj^uZ, 

where the coefficients a^ are defined in (|2.6p . 
The main result of this section is the following 

Theorem 3.1. Let us assume that z is a classical solution of N^^^z — 0. 

(i) Let us assume that a g]0, 1[, p > 10/3 and for every K <Z<Z SI there exists a constant C such that 

l|w||ci'-(K) + \\d2z\\LP{K) + \\d2XuZ\\L2(ii) + 1 1 V^2:| |i2(^) < C. 

Then for any compact set Ki CC K , there exists a constant Ci only dependent on K , C, and on the constant 
in \2.2i such that 

(ii) If in addition to the previous conditions, there exists a constant C such that 

||52X„z|U4(;,)<C, 
with a > 1/4, then for every p> 1 there exists a constant Ci only dependent on C and C and p such that 
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3.1. Lifting and freezing. The operator N^ is an elliptic (Riemannian) approximation of the sub-Ricmannian 
mean curvature operator in the right-hand side of (|1.2|) . Its linearization N^^u can be interpreted an uniformly 
elliptic operator, with least eigenvalue depending on e. It is well known that this approximating operator 
has a fundamental solution, but its estimates strongly depend on e. In order to obtain estimates uniform in 
£ we further approximate it with an Hormander type operator, a sum of squares of vector fields, which has a 
similar behavior in the direction Xi , but for which the direction 82 is the direction of one of the commutators 
( a step-three commutator!). The idea is to use a new version of the famous Rothschild and Stein lifting 
theorem, only partially inspired to the procedure in [38j . In order to deal with the non-smoothness of u, we 
will also operate a freezing: roughly speaking we approximate the coefficients of the vector field Xi with 
their first order Taylor polynomials. 

Lifting. The first step is to lift the vector fields to a higher dimensional space through the introduction 
of a new variable s. The points of the extension space will be denoted (x, s) G ft x (—1, 1) C i?^, with 
X = (xi,a;2). The lifted vector fields are defined as follows 

(3.3) Xi = d.,, + (uix) + s^)d^„ X2 = ea,,, and X3 - 9,. 
The C^'" distribution {XijX^} is a step 3, bracket generating distribution since 

(3.4) [Xi,X3] - -2sd^, and [X3, [^1,^3]] = -29,,. 

The associated homogeneous dimension is Q = 5. If a;o = ((2^0)1; (2^0)2) G $7 is a fixed point, then for all 
a; € J7, s € K one can define exponential coordinates (6^1(2;, s), 62(2;, s), 63(2;, s)), based at (xqjO), via the 
formula 

(3.5) (a;,s) = exp(^^^o)(ei (a;, s)Xi + 62(2;, 3)^2 -f 63(2;, 5)^3). 

Here for any Lipschitz vector field Z in R^ we denote by expf^^^Q^{Z) the point 7(1) where 7 is a curve such 
that 7(0) = (xo,0) and j'{s) = Z{'-f{s)). The exponential coordinates can be explicitly computed yielding 

/•I s2 ^ 

ei{x,s) = (2; - 2;o)i, ee2{x,s) = {x - ^0)2 - {x - xo)i( / u{"f{T))dT - — ), and 63(2;, s) = s. 

Jo •J 



Note that if dE{x,xo) = \/{x — xq)\ -(- (x — Xi^\ is the Euchdean distance in fl then for x sufficiently close 
to Xq and for a certain constant C > (both depending on the C^'" norm of u) one has that 

C"^d£;(x,xo) < v/ei(x,0)2-|-£2e2(2;,0)2 < CdB(x,xo). 
Next, we define an analogue of the first order Taylor polynomial of u as 

(3.6) Plu [x) =m(xo) + ei(x)XiM(xo,0) + e2(x)X2u(a;o,0), 
where ei(x) = ei(x,0) and 

(3.7) ee2(x) = (x - xo)2 - (x - xo)iu(xo). 
We remark explicitly that 

(3.8) |e(e2(x,0) - e2(x))X2w(xo,0)| = |((x - xo)i( / \u(j(t)) - m(xo)] dr)a2u(xo)| 



REGULARITY OF NON-CHARACTERISTIC MINIMAL GRAPHS IN H^ 11 

Freezing. At this point wc introduce an appropriate freezing of the vector fields by defining 

(3.9) Xi^^g = dxi + {P^^^u{x) + s^)dx^, X2^XQ = edx.^ and Xz,xo = ds- 

Observe that \X\^^^^Xj,^x^ is a distribution of smooth vector fields satisfying Hormander's finite rank 
hypothesis with step three. We denote by dx^^--,) the corresponding Carnot-Caratheodory distance and 
remark that the homogeneous dimension of the space is 5. We also need the Riemannian distance function 
'^a:o,e('j ') defined as the control distance associated to \X\_xai^2.xaT^'i,xf^- Define exponential coordinates 
(ei, 62, 63) in a neighborhood of xq through the formula 

(x,s) =exp(^^o)(ei(a;,s)Xi,^o + 62(0;, s)X2,a;o + 63(2;, s)X3,2,J. 

Note that 

ei(a;,0) = (x- 0:0)1, e.2{x,s) = -[{x ~ xq)2~ (x ~ xq)x{ \ Pl^,u{l{Ty)dT - —)\ and e3(x, 0) = 0. 

It is well known (see for instance the discussion in [71 Section 2.4] that (M^, dxo.e) converges in the Gromov- 
Hausdorff sense to (M'^, d^o)- In particular one has that for each fixed x and s, then dx^^^^^ix, s), (xq, 0)) -^ 
dxo{{x,s), (a;o,0)) as £ ^ 0. Moreover the volume of the balls i?e((xo, 0), i?) in the dxg.e metric converges 
to the volume of the limit Carnot-Caratheodory balls, i.e. \Bg{{xo,0),R)\ —^ \Bo{{xo,0),R)\ as £ — > 0. In 
particular, for £0 > £ > sufficiently small there exists a constant C > depending only on £0 such that 

(3.10) \B,{{xo.O),R)\>CR^. 

All this can also be seen explicitly in our special setting by observing that there exists a constant C > 
such that for x near xq, one has (see [32J) 

C~^dx„,e{{x, s), (xo, 0)) < y ef (x, s) + min(e|(a;, s), (£62(2;, s))^/^ ) + e§(a;, s) < Cdxa,e{{x., s), (a;o, 0)). 

/ xi/6 

C-id^„((x,s),(xo,0)) < {el{x,s) + {ee2f{x,s)+el{x,s)\ < Cdx„{{x,s),{xo,Q)). 

Recall also that for a; and xq sufficiently close, £0 > £ > sufficiently small there exist positive constants 
Ci, C2 depending only on £0 and on the C^'" norm of u such that 

(3.11) (iB((x, s),(xo,0)) < Cida;o,e((a;,s),(a;o,0)) < C2dxo{[x,s),{xo,Q)). 

Since |e2(x, s)\ ^ 00 as£ — > one has that for a fixed (x, s) thenlim£_»o dxQ.e{{x, s), (xo,0)) ~ dxo{{x, s), (xo, 0)). 
In the following we will denote by dxo,eix,xo) the quantity dxo,e{{x,0), (xo,0)). 

Lemma 3.2. If u Cz C^'"(r2) there exists constant £0, C > and a neighborhood U of xq depending only on 
the C^°' norm of u such that for all x in a sufficiently small neighborhood of xq, and for all £0 > £ > 0, 

(3.12) |^(x) - P^^u{x)\ < C4+^(x,xo). 
Proof. Fix the points x and xq and define the C^'^ planar curve 



7(t) =exp^^^ it{ei{x,0)Xi^xo + e2(x, 0)^2^x0) j : 



so that 7(1) = X and 7(0) — xq. From the mean-value theorem, for all t e (0, 1) we can find t e (0, t) such 
that 

u(7(i)) - u{xo) = ei{x,0)Xi,xoh(t)) + e2{x,0)X2,xo{l(t))- 
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Hence for all t E (0, 1) one has 

"(7(t)) - ^x>(7(i)) = u{j{t)) - u{xo) = ei(7(t))Xiu(a;o,0) - 62(7^)^2^(^:0,0) 

= ^ I ei{x,0)Xi^xgu{'-f(t)) ~ ei{'-f{t))Xiu{xo,0) 



(3.13) 

Next, note that 



(a; - xo)i{u{x) - P^^u{x))dx2u{x) + e2{x,0){X2,XQu{j{t)) ~ X2^x„u{xo, 0)) 
+e2(a;,0)(X2,^oM(a::o,0)-X2,:EoU(xo,0)) + (62(2;, 0) - e2(7(i)))^2M(a;o,0). 



|e2(a:,0)-e2(7(i))l 



(a; - xo)i / Pl^u{-i{T))dT - y + {l{t) - a;o)iw(a;o) 



< 



{x-xa)i I (P,X7(r)) - u(7(T))dr + 

1 

+ - 

e 



(3.14) 



<|4+"(x,xo)+^ 



((7(t) - xa)iu{xa) -{x- xo)iu{'y{T)))dT 
1 



{x - xo)i / (PxX7(t)) - uijir)) )di 



From the latter, from (|3.13p and from observing that X2,xo — '''^2.2;o, ^nd 

\e2[x, Q){X2,xoU{l(t)) ~ X2,xoU{xo,Q)) \ < Cd'^{x, xo)\ee2{x, 0)| 
we obtain 

(3.15) \u{^{t))-P^^u{jm 

< C{d]+"{x,xo) + \ee2ix,0)\d1:{x,xo) + \{x-xo)i\ f |Pi„«(7(r)) - «(7(r))|dT 

Jo 

For X sufRciently close to xq we have 

|ee(a::, 0)1 < Cdxo,eix,xo) and \{x - xo)|i < -. 
Integrating p.lSp in the t variable from to 1 and using the latter we obtain 

|P,X7(t)) - uijiT))\dT < C4+?(x,xo). 
Substituting this estimate back in (|3.15p we conclude the proof. 



D 



The frozen operators. The freezing process described earlier allows to introduce 'frozen' sub-Laplacians 
operators N^.xo formally defined as N^, but in terms of the smooth vector fields Xf^ instead of the original 
non-smooth vector fields Xi. Consider the operators 



(3.16) 



Ns,xoZ= ^ ai.j{VeU{xQ))Xi^xoXj,xoZ, 



where a^ are defined in p.6p . and 
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(3.17) N^^z= ^ aij{Veu{x(i))Xi^x„Xj^^^z, 

NxQ is an uniformly subelliptic operator with C°° coefficients, and N^^^o can be considered as its elliptic 
regularization, with coefficients dependent on e. The linear theory yields that both Nx^ and N^^xo have 
fundamental solutions T^a and T^^^a respectively (see [21], [32] and [4J). Since both T^a and Tf^^xoiC^O 
depend on many variables, the notation 

shall denote the Xi^^jo-derivative of Ti;^xo{C:0 with respect to the variable (, evaluated at the point Ci- 

Precise estimates for the fundamental solution Txa have been established in [32| and [1], while in [T5] it 
is proved that the fundamental solution T^^xo of N^^xq locally satisfies the same estimates as the limit kernel 
Txg, with choice of constants independent of e. These results can be summarized as follows; 

Theorem 3.3. f[l5\) Let zq G fl. For every compact set K (Z fl x (—1, 1), for every k G N and for every 
multi-index I = (zi, ...Tik) with ij G {1,2,3}, there exist two positive constants C,Cp independent of e, such 
that 

'"«' |vL.«)r.,„(-,OI<c. ,^_™^,„ . 

for every S.,C ^ ^ with ^ 7^ C; where B^{^,r) denotes the ball with center ^ and radius r in the distance 

dx„,e, fl'^'^^ex =^ -Xii.xo---Xi„,xo denotes derivatives of order \I\ = k along the frozen vector fields Xi^xo ■ V 
k = one intends that no derivative are applied on T^^xo ■ 



Remark 3.4. Let u be as in Lemma \3.2[ Notice that if we set 

Ti = Xi,u - Xi,j;„ = {u{x) - Px„u{x) - s'^)dx2, 

then since dx2 ho,s order three, and in view of Lemma \3.'A the operator TZ has order 2 — a at the point xq 
(in the sense of |23j J. Consequently, the estimates in (j3.18p do not continue to hold if the derivatives along 
the frozen vector fields V^^ , |/| = 2, are substituted by derivatives along Xf^ evaluated at the point ^ = xq. 

Since the measure of the Ball is doubling, with doubling constant independent of e, then the space 
(fix (—1, \),dxa,e) dx) is a space of homogenous type and the the following version of the fractional integration 
theorem holds, (see for instance "6^) 

Proposition 3.5. If K is a fixed compact set K izci Vt x (—1, 1), if f G L'^{Vt x (—1, 1)) supported in K, 
and for each (a;, s) G fi x (— 1, 1) we set 

r dP Mx,s),(C,(t)) 

J \Be{{x,s),dxo,s{{x,s),{C,a))\ 
then there exists a constant C , depending on K, r, q but independent of e such that, 

II^p(/)I|l-(ox(-i,i)) < C\\f\\L<,(ny.{-i.i)), 
and where 5 — po > and r = , '^ . 

Corollary 3.6. Let f e C^{^) and extend it to a function on Vl x (—1, 1) by setting f(x, s) := f{x). Let 
/C : (il X (—1, 1))^ -^ M. be a kernel satisfying 



\mx,s),{Ca))\<C- 



<,,((x,g),(C,a)) 
Be{{x, s), dxo,e{{x, s), {C, a))\' 
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Set u{x, s) to be defined as 



u{x,s) = / IC{{x, s), {C, (T))f{^)dCda. 
If u{x, s) — u{x) for all values of s then 

\\n\Win)<C\\f\\L,(n). 

with r and q as in the previous proposition. 

As a consequence of the definition of fundamental solution one has the following representation formula: 



Proposition 3.7. Let u be a fixed smooth function defined in ft, let N^^^ be as in p.2p . and let z be a 
classical solution of N^,uZ = .g G C°°{^). Extend both u and z to be functions defined on 51 x (—1,1) by 
letting them be constant along the s variable. For any (f> G C^{Cl x (—1,1)), ^ G 57 and s G ( — 1,1), the 
product z{^)4>{^, s) can be represented as 



2:(C)0(C,s) = / r((C,s),(C,cr)) zN^^g(t)+ ^ a,j{xo)(^X,^,_„zXj^,_g(j) + Xj^,^„zX,^^„(t)j dC 
flx(-l,l) \ *^"^ 

r,,,„((e, s), (C, a)) g(C) </>(C, fT)dC'^a+ 



Ox(-l,l) 
2 



(3.20) 



- X! / ^e,xo{{C,s),{C,(^))[a,j{xo) -a.,j{C)jX,^uXj,uz{C)(l){C)dCda 

'^ = ^Ox(-l,l) 



2 „ 

E°y(^o) J r,,,„((f,s),(c,a))(7.(C)-P,XC)-^')52X,-„z(C)0(C,'T)dCda 

■'^-^ Ox(-l,l) 

2 „ 

E«»i(^o) / ^^.=»o(C,5)r,,,„(., (C,cT))(«(C) - P^AO - a^)d2z{04>{C„(T)dCda 

7 — 1 

Ox(-l,l) 
2 „ 

E °»i(^o) j r,,,„ ((e, s), (C, fi)) (^(C) - <"(C) - ^') ^2^(0 x,,,,M. ^)dCda. 



'-^ Ox(-l,l) 



/n order to simplify notations we have set 

ay(C) =ay(V£u(C)). 
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Proof. In view of the definition of fundamental solution, we have 
(3.21) 

z{0^{i,s)^ J r,,,„((^,s),(C,a))7V,,„(z<^)(C,a)dCda 

f2x(-l,l) 

r£,xo((^, s),(C,cr)) zN,x,-,(j)+ ^aij{xo)[Xi^xgUXj^x,,(t> + Xj^xgZXi^^g<l)j d^da 
Ox (-1,1) V *^ = ^ / 

re,.o((^, s), (C, a))N,uz{C) </)(C, a)dCda 
nx(-i,i) 

+ J r,,^M,s)AC,^))iNe.o - Ne^.)z{C,<J)HC,<J)dCda. 

nx(-i,i) 

Let us now compute the difference between the operator N^u and its frozen counterpart. We will emphasize 
the presence of the variable a in those terms where the coefhcients of the vector fields involve that variable. 

2 



2 

ij = l 
2 

ij=l 
2 

- 2_^ 0'ij{xo)UXi^u — Xi^xo)Xj.u + Xi^xo{Xj,u — Xj,xo)jz{(, a) 

2 

(3-22) = ^ [a.,j{xo) ~ dijiOjXi^uXj^uziC)- 

ij=l 
2 

- Y. d,,{xo)[5.a{u{0 - P,XC) - a')d2X,.^^+X,,xo{S,Au{C) - P^,u{C) ~ <j')d2))z{0 

2 

ij=l 
2 

1=1 

2 

- J2 afi(^o)X.,.o ((^(C) - P^MO " <^')d2)z{C). 

The integral of the last term in (|3.22p becomes 
2 „ 

(3.23) E«^i(^o) j r,^xo{{ts),{C,a))X,^xo{{u{0-P^,u{0-a')d2)z{C)^{C,cT)dCda 



Ox(-l,l) 
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(integrating by part) 



2 r 



(3.24) 



Ox(-l,l) 



E ««i(^o) / re,.o((e, s), (C, 'j)) HO - P,XC) - <j')d2z{0 ^.,.o</'(C, '^)rfCrf'T 



nx(-l,l) 



Inserting all terms in the preceeding formula we conclude the proof. 



D 



Lemma 3.8. Let u and xq Cz Q be as above. There exists a neighborhood U of [xq, 0), possibly depending on 
e, such that for all (^, s) Cz U one has 



(3.25) 



Xk,xoXl,xoii,s) - Xk^xoXl,xoixo,0) ]^xo,e^E,xoi-, (C; cr)) 



<d,M,s),i^O,0)){\yio,e^e,xo{{i,s),iC,^))\ + \^Le^e,xoiixo,0),iC,a))\], With \I\ ^ 4. 



Proof The proof follows from the mean value principle: Set L — dx„{{^, s), (xq, 0)). For every (5 > consider 
a horizontal curve 7 : [0, L + (5] ^ fi x (— 1, 1), parametrized by arc-length and joining (^, s) to (xq, 0). For 
every f E C^ one has 



\f{ts)-f{xo,0)\< J^ \-f{^{t))\dt^ j \{Vf{-i{t))n'{mdt<{L + 5)sM^^o,ef\. 



There exists neighborhood (xo,0) G U depending on the C^ norm of / for which we have |Va;„_£/(C, cr)| < 
2(|Va;j,^e/(a:;o, 0)1 + \V xa.ef{^is)\) for all (C, cr) G U. The lemma now follows from choosing 



/ ^ Xk^xQXi^xaXi,xcXe,xo(,'i iC^'^)) 



and observing that the smoothness of T^^xo depends on e > 0. 



D 
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Lemma 3.9. Let z be a smooth solution of N^,uZ = 0. For any s G (—1, 1) and xq G fl, one can represent 
the 'frozen' second order horizontal derivatives of z at xq as 

(3.26) 

Xk,xaXl,xa{z{xQ)(l){Xf), s)) = 



nx(-i,i) \ *^=^ 

2 „ 

+ X! / Xk,xaXi,xo{xo,s)T^.xo{-AQ^'^))\^ij{^^)^0'i]{0)Xi,uXj^uZ{C,)(t){(^,(j)d(^da 

2 „ 

-^ayK) j Xk^x,Xux,X^^,s)T,^x,X-,{Q,(^))[u{C)-Plu{0~a^y2Xj.uz{C)4>{C,(T)dCda 

^=^ nx(-l,l) 

2 „ 

+ ^aii{xo) / Xk^xoXl^xoX^,xo{xo,s)T^xoi■,{C^f^)){u{C) ~ Pxo'^iC) ~ <^'^jd2z{C) (l){C„a)dCda 
'=^ nx(-i,i) 

2 „ 

+ ^ aa(a;o) / Xk,xaXi^xo{xo, s)Te,xa{-, (C, cr))(u(C) - Pxo^iC) - cr^)d2z{C) X^^xo<P{C,, s)dC,d(j 



'-^ Ox(-l,l) 

Proof. Since the proof is similar to that of [171 Proposition 3.9], we only sketch the argument for the most 

singular term in the representation formula, i.e. 

(3.27) 

2 „ 

/((^,s),(a:o,0))=^a,;i(a:o) j X,,,„(e, s)r,,„(., (C,a))(u(C) - P^XC) - ^2)92^(0 </)(C, s)dCd'T. 

*=-^ nx(-l,l) 

We want to show that 

(3.28) Xfc,,„Xi,,„(a:o,0)/(-,(.To,0)) = J^^), 

where 
(3.29) 

2 „ 

/(2)(a;o) =^aa(a:o) j Xfc,,„X,.,„(a;o, s)X,,,„r,,„(., (C, a))(w(C) - P,>(C) - a2)a2^(C) 0(C, s)rfCd^ 
*=^ nx(-i,i) 

Note that the latter is well defined in view of the estimates p.l2p and p.lSp . To show (|3.28p we con- 
sider a family of smooth test functions Xxo,e{{£,Ts),{C,,a)), satisfying for some choice of C > and for 
small e > 0, (i) < x.o,e((e, s), (C,^)) < 1,' (ii) X.o.e(((C,s), (C,'^)) = if d,„,e((C,s), (C,a)) < 2Ce, (iii) 
Xxo,e((e,s),(C,a)) = 1 if d,„,e((e,s),(C,a)) > AGs, (iv) |V,%„Xxo,e(((C, s), (C^))l < Ce-\'\ for aU multi- 
indeces /, for some choice of C > and for small e > 0. For the existence of such function see [T^] and note 
that the construction argument in that paper uses only the estimates on the fundamental solution. Define 
the smooth approximation 

2 r 

(3.30) /,((e,s),(xo,0)) = ^au(xo) j X,,,„(C, s)r,,,(., (C, a)) 

'=^ Ox(-14) 



(u{0 - PluiO - a^) X.O.M, s), (C, ^)) d2z{0 </)(C, s)dCdc 
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Once we establish the bounds 



(3.31) 



/,((e,5),(xo,0))-/((e,5),(xo,0)) 



<Ce", 



foiQ (ij:o_e((^, s), (xo,0) < e, then invoking [T71 Proposition 3.2] we immediately conclude (|3.28p . To prove 
the first estimate in p.3ip we first observe that for e > suSiciently small 



(3.32) 



4((C,5),(xo,0))-/((C,5),(a:o,0)) 
<C 



•d.o.a(C,<^),(a;o,0))^+"|92^(C) </'((, 5)MCda. 
Since 

d.o,e((C, ^), (a^o, 0)) < C(d,„,e((C, ^), (e, s)) + d,„,,((xo, 0), (e, s)) < C(d,„,e((C, ^), (C, s)) + e), 
then we conclude 

<,,((?, s),(C,-7))2+" 



4((C,s),(xo,0))-/((e,s),(xo,0)) 



<C 



Next we turn to the second estimate in p.3ip . Observe that 

Xk,.oXi^.o{OlM. s), (xo, 0)) - /(2)(^o) 



<Ce 



2+a 



where 



(3.33) ^1= J 



ox(-ia) 



(3.34) ^2= I 



Ox(-14) 



<|Ai| + |A2| + |^3|, 
^fc.a:o-'^/,a:o(?, «) " Xk,x„Xl^xo{xQ, s) \ Va;o,er(-, ((, a))- 

■ (uio - p^AO - ^') x.oAi^, 5), (C, ^)) a2z(C) 0(C, ^) 

^fe,2;o^i,2;o(a;o,s)Va;o,er(-,(C,cr))- 

• (uio - P.XC) - ^') (i - x.o,e((C, s), (c, ct))) a2z(C) </-((, '^) 



d(^da, 



d(^da, 



Hn view of the dependence of the neighborhood U from e in Lemma 13.81 possibly one needs to work with a smaller scale 
da;o,e((5, s), (a:o,0) < o(e) depending on the C^ norm of ^s,xq 
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and 



(3.35) ^3 



r2x(~i,i) 



^i,xoV,„,,(C,s)r(.,(C,a))- 



(m(C) - P,XC) - ^2) x,,,„(e)x.o,s((e,s),(C,'T)) a2z(C)0(C,'T) 



dC,d(j 



Slx(-l,l) 



V.o,s(^,s)r(.,(C,a))- 



(w(C) - P,XC) - ^2) Xz,,„X,,,„(e,s)x.o,e(-,(C,^)) d2z{Q)cj>{C,a) 



d(^da 



Invoking Leinnia l3.8l one can complete the proof arguing as in [171 page 734]. As usual, we examine in detail 
only the integral Ai which contains the most singular integrand. Note that if C, Gsuppx^.^ g((^, s), •) then 
c^a;o,e((?,s),(C,cr)) > 2Ce, consequently 

dxo((2;o,0),(C,CT)) > ^rfa;o,e((^,s),(C,CT)) - d^„^e{{^, s), {x^,0)) > E. 

In view of Lemma I3.8i (13.31) and Lemma 13.21 we have 



dQda 



(3.36) \A,\<C f d,oms),ixoMdit"ii^o,0)AC,^))HO ,^..'^^"!^J^^;^'^^^ ,,,, 

id,o((:Eo,o),(c,a))>e |S((a;o,0),rf:ro((a;o,0),(C,cr)))| 



C I d^M: s), {xo,0))dl+"{{xo, 0), (C, a))0(C) ,r,..fi^^:%f^X\,,, dCda < Ce". 

'd^„ai,s),{c,-^))>e l^((C,s)>c?xo((C,s),(C,a)))| 



n 



Using the representation formula above, the fractional integration result in Proposition 13. 51 and Corollary 
3.61 we finally can proceed to the proof of the main result of the section: 



Proof of Theorem \3.1\ We will prove (i) only. The proof of (m) follows along a similar argument. Using the 
representation formula p.26p one can represent the second horizontal derivatives Xk^xo^i,xo{z{xo)(j){xQ, s) of 
z at any point (xq, s) e 57 x (—1, 1) through integral operator with kernels of the form 



X 



k^x„Xi^x„X,^x^{xo, s)T^xoi-, (C, cr)){u{C) - Pl-,u{C) - cr^j 
Xk,xoXi^xo{xo,s)r^^xo{-, {C,<^))[dij{xo) - ay(C)j 



and 



Xk^xoXLxoixo,-s)r,,x„{;{C,<T)){uiO - P^AO ~ '^^) 



To establish the non-singular character of such kernels one needs to invoke the estimates on the derivatives 
of the fundamental solution of the frozen operator T^xo in Theorem 13.31 To estimate the L^"/'^ norm of each 
term in the right-hand side of (|3.26p one uses the fractional integral estimates in Corollarv l3.6l The 'worst' 
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possible term is the one corresponding to three derivatives on T, i.e. 

2 . 

(3.37) ^aa(xo) J Xfc,,„X,.,„X,.,„(xo, s)r,,„(-, (C, <t))(u(C) - P,>(C) - a^^aa^CC) 0(CKd<T 

'=^ Ox(-l,l) 

= JlC{{xo,s),{C,a))d2z{(:)HC)dCda, 
with (in view of Remark 13. 4p 

/C (xo, s), (C, cr))\ < C—— — \- J- s .. .,, ■ 

|Se((a;,s),d^„,e((xo,s),(C,o-))| 

Note that the expression in p.37p does not depend on s. Moreover since from the assumptions one has 

p > 10/3 >, and hence p > 50/(15 + 10a) then CoroUarv 13.61 vields immediately that the integral in (|3.37p 

is in L^ with 

5p 10 



5 — pa 3 
The rest of the terms in the right-hand side of (|3.26p are estimated similarly. 

At this point we have proved that the function xq -^ Xk^xo^i.xozixa) is in L~{K). Next we observe that 

{0 if ; = 2 and fc = 1,2 

u{xo)dx2u{xfi)dx^z{xQ) if ; = 1 and fc = 1 
if ; = 1 and fc = 2 

In view of the hypothesis dx2Z e LP{K) then one finally concludes 

Xk^uXi^uZ - Xk,x„Xi,,,{; 0)z e LP{K) 

withp> 10/3. D 

4. CACCIOPOLI type INEQUALITIES: WJ"l'oc A PRIORI ESTIMATES. 

In this section we prove a a priori estimate, in the Sobolev spaces W^™o^(^) for solutions z of equation 
(|2.8p . under the assumption that u is smooth and satisfies (12. 2p . 

The starting point is a Caccioppoli type inequality for derivatives of solution of (|2.8p in the directions Xi. 
By Lemmas 12.61 and 12.71 these derivatives solve the same equation (|2.8p . with a different second member, 
hence we will focus on this PDE. 

Lemma 4.1. Assume that /o € Lj^^{fl), and that z g W^'ig^{fl) D W'l ;oc(^) *'' "^ solution of equation Ii2.8\) 
For every p > i there exist constants Ci — Ci (p, M) with the constant M as in i2.^) and independent of e 
and z such that for every non-negative (j) € C^(f2), we have 

|V,(|.|(P-1)/2)| V < c, {^j |zrl(</.2 ^ |^^^|2) _ j /|,|P-3^^2^ 

Proof. Let us multiply both members of equation (|2.8p by \zY'~^z(fP' , and integrate. We obtain 

(4.1) j f\zr'zc^^ = j xA^M^^X,z)\zr^zc^^ = 

J J ^^Jl + \VeuY ' 

(since XI = -Xi - diu, X^ ^ -X^) 

= -[d.u6a ^'^^"'f x.zlzr^zc^^- 
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HP - 2) / -^ig^X.zX^zr^^ - 2 f -^ig^X,z\zrh4'X.^. 
J Vl + I^e"r J vl + IVewr 

This obviously implies that there exists a constant C > such that 

(4.2) fc^ J |V,(|z|(^-^)/^)P<^^ < C| |zr i(|V,</-p +<P')-J f\zr'z<^\ 

concluding the proof. D 

Lemma 4.2. Let p > 3 be fixed and u be a function satisfying the bound (12. 1|) . Let f G C°°{^), and let z be 
a smooth solution of equation i2.8\) . There exist two constants C and C which depend on p and the constant 
M in i2.<^) but are independent of e and z such for every (j) G C^(0), </> > 0, 



(4.3) / |Ve(|VeZ|(P-l)/')| V'P <C (|Ve0|' + (l>Y + 



Proof. Since z is a solution of equation (|2.8p then by Lemma 12.61 si — Xiz satisfies equation 

M,s, = X,{f) + xj ^^^^^^d2uX2z] - 

_xJxJ^ME^^)x,z]+d2uxJ^MB^X,z] 
Using Lemma |4. II we deduce 

(4.4) / |Ve(|si|(P-l)/2)|202p < ^^ y |^^|p^l^|^^^|2 ^ <^2)^2p-2_ 

aj2(VeW) 



/i„r=,.x.(x,(-il 



d2uX2Z 6^P~ 






a,,(VeW) ^^.^),2p 



- /|.ir-^.i92^xj '!''^^f\ x,z]^'P~ f\s^r's,X,fc^^P = 
(integrating by part all terms in the right hand side) 

= Ci j \s,r\\^A\' +<^')<^^P-^ + I \s,r's,^^^^{d2ufX2Z<^'P + 
+ (p-2) [ \s,r^X,S,^S^^d2UX2Z<f'P + 
+2p f\s,r's, ''''^^'''\d2uX2Z<f>'P~'X,cl^+ 
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IP-3, a „, a2j(Ve 






(5 „ 

l|PrA2P-^ / IV.r|P^2p^ 



+ C / \X2{d2u)\P^'P + I \V,z\Pcj,^P + ^ I |,SirV202p^ 



+ ^ /■ |J|(2p+l)/7^2p ^ f |^^|(2p+l)/5(|y^^| ^ ^)(2p+l)/5^(8p-l)/5 < 

< Ci I (I V,0P + 02)'' + 11 I V,zr,/.2p ^ ^ I |,^|P+l/2^2p^ 

+ J \X2{d2u)\Pcj)^P + (^y |sir3|V,Si|202p^ 

iv^«iiVezrv'^ + y I iv^^nvezrv'^ + J iv^^iivezr v'^~'ive0i+ 

+ £ /■|/|(2p+l)/7^2p^ /"|/|(2P+l)/5(|V,(/,|+(/,)(2P+l)/5^(8p-l)/5 

It follows that 

(4.5) I |Ve(|.i|(^-i)/2)p02p < c(| (iVe^p + 0^)^+ 

+ / |VeZr+l/2^2p ^ f \X2{d2u)\Pcj)^P + I |/|(2P+l)/7^2p ^ /^ |^^| (2p+l)/5 (|y^^| ^ ^) (2p+l)/5^(8p-l)/5 

An analogous estimate holds for S2 ~ X2Z, i.e. 

(4.6) / |Ve(|.2|(^-^)/^)p0^'' <c(f (|Ve0p + 4>'Y+ 
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|V?u||V,z|P-U2p+ /"|V?uP|V.z|P-V'P+ /"|V?u||V.z|P-V'P-MV, 



The conclusion follows immediately. 

D 

Theorem 4.3. Let p > 3 be fixed and u be a function satisfying the bound (12. ip . Let f € C°°{fl), and let z 
be a smooth solution of equation 112. 8\) . There exists a constant C, which depends on p and the constant M 
in 112. 2fl but is independent of e and z such that for every (f> G C^{fl), </» > 0, 

I |V,zr+l/202p ^ J |^^^|^^^|(p-l)/2) |2^2p < 

< C[J\Z\^P+^^^P + j {<t>^ + \^e<t>?Y + j \Ve{d2u)\P(p^P + 
+ / ('|/|(2P+l)/702^ + |/|(2P+1)/5(|V,0| + 0)(2p+l)/5^(8p-l)/5\ ^ j |y2^|2p^2p\ 

Proof. By Proposition l2.2[ and Lemma W?I\ calling si = Xiz we have 

f |sir+l/2^2p < ^ /■ (^|^|4p+2^2p ^ |^|(2p+l)/2^2p ^ |^| (2p+l)/2 |y ^^| (2p+l)/2^ ^ 

+<5 f {\V,z\P+'/^ + \X2id2u)\P)<l>^P + fi\Vlu\^P + \Vlu\P)<P^P+ 



+C / (|/|(2?'+l)/VP + |/|(2P+l)/5(|V,0|+0)(2P+l)/50(8P-l)A 



To estimate S2 — X2Z we argue in the same way and obtain 

C 



|52|P+V2^2p< ^ / (|z|4P+VP+|z|(2P+l)/202p + |^|(2p+l)/2|y^^|(2p+l)/2^^ 

+5 /(|V,z|P+i/2 + |X2(a27/)r)02P+ fi\V^,u\^P + \Vlu\P)<l)^P+ 



+C f (|/p?'+l)/702p ^ |^|(2p+l)/5(|y^^| ^ ^)(2p+l)/5^(8p-l)/5j 

Hence, if 6 is sufficiently small 

C 



(|si|P+l/2 + |52r+'/')<^'^ < y / (|z|^P+'02P + |z|(2P+l)/202p + |^|(2p+l)/2|y^^|(2p+l)/S 

+ J (I Ve0p + 0^)" + J \Ve{d2u)\Pcj,'P + j{\Vlu\^P + \Vlu\P)cf'P + 
+C I (|/|(2P+l)/702p ^ |^|(2p+l)/5(|y^^| ^ ^)(2p+l)/5^(8p-l)/5^ ^ 



The conclusion follows from the latter, (|4.5p and (14. 6p and the Holder inequality 

|V,uP(i2P< /■ |V,M|2P,i2p ^ /■ ,2p^ 



D 
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Next we iterate once the previous result 

Theorem 4.4. Let p > 3 be fixed and u be a function satisfying the bound (|2.ip . Consider a function 
f G C°°{Q), and z a smooth solution of equation 112. 8\) . Let fli,fl2 so that fli CC 5^2 CC fl. There exists a 
constant C, which depends on p, on fli, and on the constant M in i2.S\} . but is independent of e or z such 
that 

II ||P+l/2 , >r^ III / |(p_i)/2||2 < 

|/|=2 

Here L is a multi-index and V^z denotes derivatives of order \L\ along Xi_^, i — 1,2. 

Proof. If z is a solution of p.8|) then, by Lemma [2.71 the function S2 = X2Z is a solution of the equation: 

M,S2 = fo 
where 

h=X2f-xJ «'^(^^") 



-X, I X2 ( "^■'"^^^''\. ) X,z I - d2uX2 I -;M^^^^^X,z 




Let us choose il^ such that ili CC ils CC il2- By Theorem 14. 31 there exists a constant C independent of e 
such that 

(4.7) ll^^li;;'^?V2(,3yll|Ve.2|(-^)/^||^.2(,3,< 

< c(ii/oiii(C.\v;.(o.) + ii-ic-^(o.) + ii-ii';.^-(..) + ii^2ii'/--(02) + 1 

We note that for three fixed functions /, g, h, 

ll/5'l|lL(2p+l)/5(n2) - 

<-||f||2p/3 II ||4p(2p+l)/4p-3 , ||. ||4p(2p+l)/4p-3 

^ \\J IIl2p/3(Q2) "T I|y|li4p(2p+l)/4p-3(s^2) ^ I I ""I I L<lP(2p+l)/<lp-3 (Slj) " 



< ii/ii'^,/^.. . + mttLr.^. + \\h\\%tl.^_. + C. 



If follows that 



L2p/3(n2) + II5|Il4p + 2(02) + II'1|Il1p + 2(S12) ^ 

mT ||(2p+l)/5 

||/0||j^(2p + l)/5(Q2) - 

||.||(2p+l)/5 II 2 y ||(2p+l)/5 II ||(2p+l)/5 

^ \\J ll^l,(2p + l)/5(J^2) ^ II ^£"'^^e^llL<2P+l)/5(n2) ll^e^^e^llL(2p + l)/5(Q2) 

+ IIV3mV zll^^P+^)/^ +llV2uV2zl|(^P+^)/^ < 

T|| VeUVeZ||^(2p + i)/5jj^^j + II ^£"^e^llL(2p + l)/5(Q2) - 

< ^A|f||(2p+l)/5 II ||4p+2 



ill ||4p+2 , II ||2p/3 II ||4p+2 , II ||2p/3 

+ ||M||^2,4p+2(jj^^ + Il"llw'i^'2''/^(n2) "^ Il^llw^^''+"(n2) ^ ll^llM/f'''/^n2) ^ 

Arguing in the same way with the function si = Xiz we conclude the proof. D 



Iterations of Theorem 14.31 yield the following 
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Theorem 4.5. Let p > 3, m > 1 be a fixed positive integer and u be a function satisfying the bound ()2.ip . 
Assume that f G C°°(f2), and let z be a smooth solution of equation 112. 8\} in Q. If ^i CC ^^2 CC fi then 
there exists a constant C which depends on p, Qi and on M in \2.'2\) . but is independent of e or z such that 
the solution satisfies the following estimate 

II ||P+l/2 , y^ my/ |(p-l)/2||2 < 

\I\=m+l 



<^(ii/iep,^c^v.(,^,+iH 



|(2p+l)/5 _LlU,l|4p+2 

i,4p + 2/ 



_L|U||4P+2 , I 



U||4p+2 +|U|pP/3 +|U||4p+2 -|-||z||2P/3 , 1 



5. A PRIORI ESTIMATES FOR THE NON-LINEAR APPROXIMATING PDE 

We now return to the equation L^u = 0. Let u be a smooth solution satisfying (|2.ip . In view of Proposition 
12. H and Theorem l3.1l (i) we have the foUowing statement: for every open set Q.i (ZCl Vt there exists a positive 
constant constant C which depends on Vti and on M in (|2.2p . but is independent of e such that 

(5-1) \\'^\\wi'^"'\a^) + II^2u|Im.^i,2(o,) + Il"llc],'°(f2i) ^ <^- 

Our first step is the higher integrabihty of the Hessian of u. The proof rests on the estimates obtainted 
from the freezing technique in Theorem 13.11 and from a new EucHdeaio CacciopoH inequahty ([H 



Lemma 5.1. Let u be a smooth solution of 

LeU = 0, 
in fl C M.^ satifying (j2.ip and denote v ~ d2U. For every open set fli CC fl, for every p > 1 there exists a 
positive constant C which depends on fli, p, and on M in h2.^) . but is independent of e such that 

Proof. In view of Lemma 12.51 the function v = d2U satisfies the equation: 

with 

Hence, applying Lemma 14.21 one has 

(5.2) I |Ve(|Ve«|(^-l)/^)p0^P <C[J (|Ve</)|^ + C^Y + 

|j|(2p+l)/5(|y^^| ^ ^)(2p+l)/5^(8p-l)/5^ 

There exist positive constants C\ — Ci(|Ve0|, 4>, M) and C2 = C2(M) such that for p = 3 we obtain 



rather than subelliptii 
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(5.3) / I V^vp^e < Ci + C2 ( /" I VeWp+l/'0'^ + /"(I + I VeV| + I V2u|)^/5^23/5(|y^^^| ^ , 



^7/5 



+ 



It follows that 

(5.4) 1 1 V?.|^/ < ^ 1 1 V?.| V^ + -5 / I V,.r/ + ^- 

Analogously, if we set z — Xiu, or z — X2U, using Lemma 12.41 and arguing as above we have 



2.,|2 i,6 ^ ^2 / |Y72„,|4 16 , ^^ , n / IY7 „,|3i6 



(5.5) y I v^zi^c/." < ^ y I v^uiv-^ + T + ^V '^^'' 

Using Lemma [^31 (15.41) and (|5.ip . we obtain immediately 

Hence 

(5.6) J\VMU'<Ci + C2j\VluU' 

Consequently, from the latter and (|5.5p we deduce that 

(5.7) yiv^ziv<ci+c2yiv2wiv 

Next, from the intrinsic Cacciopoli inequalities (j5.6p and (15. 7p we deduce an Euclidean Cacciopoli inequal- 
ity: Note that 

\VeXiz\ < \Xfz\ + C2\d2Xiz\ < \Xfz\ + C2\vd2z\ + C2\Xid2z\ < 

(since d2Z — d2Xiu — v'^ + Xiv) 

\Vlz\+C2\'^lv\+C2\VM+C2- 

From the latter and (|5.6p and (|5.7p we infer 

(5.8) J |VsV,zp06 < C2 (y IV^H"/ + y IV^zp^e + l) < C2 y |Vez|V' + C7i 

Now we can apply the standard Euclidean Sobolcv inequality in R^ and obtain 

( y (|VeZ|03)6^ 1/3 < ^^ y |^^(^^^^3)|2 < ^^ y |V,z| V _^c,< 

(using Holder inequality ) 

1/3 / /■ s 2/3 



C2{ U\Wez\cl>r) (/ |V.z|3) +Ci. 



By (|5.ip and the fact that |Ve2:| < |V^u|, we already know that |Vez| G Lf^^. In fact 

( / \^eZ\'f' < ( / \W,zr/'f'\supp{^)\'/''. 

<y supp{(l)) J supp{(p) 

Recall that C2 doen not depend on |V£(/)|. If we choose the support of ip sufficiently small, we can assume 
that the integral J^^ ,,. |Vez|^ is arbitrarily small. It follows that 

{fi\V,z\^ry'<C, 
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and consequently, by (15. 6|) 

|4j,6 



But this implies that \S/ E(y eu)\ < |V^u| + |Veu| + u^ G Lf^^. This imphcs, buy the standard Euclidean 
Sobolev Morrey inequality in K^ that 

By Theorem 13.11 (ii) it then follows that for every r > 1 there exists a constant C > independent of e 
such that 

D 



In order to boostrap regularity we apply Theorem 14.51 to the non linear equation L^u = 0, and obtain 
immediately the following: 

Lemma 5.2. Let u be a smooth solution of L^u ~ 0, satisfying p.ip . Set z — Xiu, i — 1,2 v — d2U. For 
every open set fli CC 5^2 CC fi, for every p > 3, and every integer m > 2 there exist a constant C which 
depend on p, m fli and on M in i2.2\) . but is independent of e such that the following estimates hold 

(5.9) iNi;r:av.(,^,+iHi:;:a..,,^,< 

(5-10) Mwr--^in.^ + Mwr-^-in.)< 

Proof. In view of Lemma 12.51 the function v solves an equation of the form 

MeV = f^, 

with 



5.11) ft, = =i; — 3 — =i'A,-t; — Aj — ^= 

Analogously, the function z = XiU solves the equation 

MeZ = /„ 

of the form with 

(5.12) /, - -[Xk,xM^^^^=) - xJ^ME^^[Xk,X,]u 

Hence 

(5.13) |/,| + |/,|<C(l + |Vez;|+Vez|), 

for some constant C depending only on M in (|2.2p . 
Applying Theorem 14.51 to z and to v, yields: 



V . 






|P+l/2 +|U||P+l/2 j_ Y- |||V/.|(P-1)/2||2 ,„ _^ \- |||V/„|(P-1)/2||2 



|(2p+l)/5 



- *" ("•^^"^,""''""+''^"(02) "^ ll/''llvv'--i>(2p+i)/5(s^^) 
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ill ||4p+2 I II ||4p+2 I II ||2p/3 

(substituting (|5.13p in the latter) 



<c(||v,z|iL^i+i)/L.„,„, + llv..|i(^^+^)/^ 



i-l.(2p+l)/5,„ ^ I II V ei'll „,_i,(2p+l)/5 



(using Holder inequality) 

^r'fw l|4P+2 I II ||2p/3 I II ||4p+2 _|_|| ||2p/3 

We have proved (|5.9p and substituting p = 3 yields (|5.10p 



D 



The main result of this section is the following a priory regularity estimates for solutions of the approxi- 
mating non linear equation: 

Theorem 5.3. Let u be a smooth solution of 

L^u = 0, 

in fl CZ M.^ , satisfying (|2.ip . For every open set fli CC fl, for every p > 3, and every integer to > 2 there 
exist a constant C which depends on p, m i7i and on M in ^2.2fl . hut is independent of e such that the 
following estimates holds 

(5.14) l|w|kr-(ni) + \\d2u\\wr-''{n,) < C. 
Proof. The proof follows from the estimate 

(5.15) ||Xlu||^y„.-l,p(J^^) + ||52u||^y™-i,p(j^^) + \\Xiu\\^rn,2f^^^^ + \\d2u\\^^,2f^^^^ < C, 

which we prove by induction. 

First step: m = 2 

By Lemma 15.11 we already know that there exists a constant such that for every p 

(5-16) IHC-(..)^^- 

We need to show that v £ WjJ^^ for every p, and that Xiu, v € W^l^^. 

Note that we can not yet invoke Lemma l5.2| since it only apply to higher order derivatives. 

Recall that w is a solution of M^v = f^, where /„ is defined in (|5.1ip By Theorem 14 . 31 there exist a constant 
C such that 



|V,«ri/^^2p^y |^^(|^^^|(p-l)/2)|2^2p< 
(|/,|(2P+l)/702p ^ |/„|(2P+1)/5(|V^0| + 0)(2P+l)/50(8p-l)/5) + f |v2u|2P02p^ < 

(using (J5T3)1 . (1511) and ({511)1 ) 

<C{1+ f |VeZ;|(2P+l)/5(|Ve0| + 0)(2P+l)/50(8p-l)/5 + / \\/ ^ {d2u)\P d^A < 
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< C(l + / |V,w|P(|V,0| + ^)(2P+l)/50(8p-l)/5^ . 



From Lemma |5. II the right hand side is bounded for p — 4. An obvious boostrap argument yields V^u € L^^^ 
for every p. Moreover, choosing p = 3, we also infer 



V2« G Ll, 



To conclude the first iteration step we observe that the function z ~ Xiu solves M^z = fz, where fz is 
defined in (|5.f2p . Theorem 14.31 and estimate (|5.f3p yield that there exists a constant C such that 



j\v,z\p+^'^<t>^p+j\^,{\y,z\^p-^^'^)\'<t>''p<c. 



Choosing p = 3 we obtain Xiu ^ ^eioc 

Main iteration step: to > 2 

Assume (|5.f5p holds for for a fixed value of m. 

Let 512 be as in Lemma 15.21 In view of that result we infer 



KcfUvW^P^'^ _L|U||4p+2 _L|U,||2p/3 i|U||2p/3 mA< 

by induction assumption 

i/7'P 



The same boostrap argument used above implies v, z d ^eioc ^'^^ every p. 



Invoking (fSTOll 

(5-18) Mwr-^^^in,) + Mwr-^^^in,)< 



<c{\\v\ 



\^\\wr-'-'H^2) + Il"llw'r''(n2) + ll^llwr''(02) + ij ^ C'' 



concluding the proof. D 

6. Estimates for the viscosity solution 

In this section we turn our attention to the proof of regularity for vanishing viscosity solutions u of 
equation p.2p . The regularity is expressed in terms of the intrinsic Sobolev spaces Wq '^(il) and rests on the 
a priori estimates proved in the previous section in the limit e — > 0. 

Let w be a vanishing viscosity solution, and (uj) denote its approximating sequence, as defined in Definition 
11.11 For each Ej and function Uj we set Xij = di +Ujdx2, ^2,] = £jdx2 the corresponding vector fields, and 
let Ve and W^!P{Q) denote the natural gradient and Sobolev spaces. We also let u , Xi = di + ud2, and 
Vo = (Xi,0) denote the coefficients and vector fields associated to the limit equation and the limit solution 
u, while W^'P^fl) will be the associated Sobolev space. Note that Vb and W^'P{^) are the usual gradient 
and Sobolev space. 

Theorem 6.1. Let u g Lip{VL) he a vanishing viscosity solution of (II. 2[) . and set Vj — d2Uj. For every ball 
B{R) CC ft and p > 1 there exists a constant C > such that 

(6.1) l|Ve,%lllV^'P(S(i?,)) + WVjWL^-iBiR)) + \\Vj\\w^fiB{B,)) ^ ^ 

and 

(6.2) XijMj -^ Xu, X2jUj -^ 
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as j -^ +00 weakly in W^' i^^^fl). Moreover equation U.^) can be represented as 

X^u = 
and is satisfied weakly in the Sobolev sense, and hence, pointwise a.e. in fl, i.e. 

XuX*(t) = for all G C^{Vl). 

Proof. The uniform bound on HwjUiooj-^^^j)) follows from the definition of vanishing viscosity solution. The 
bound on ||wj||^yl,2,^(.^^^ is a consequence of (j5.14p . To prove the remaining estimate observe that for any 

function w: 82X1 jW = Xijd2W + d2Ujd2W. Substituting w — Uj and in view of (|5.14p we see that there 
exists positive constants Ci,C2 depending only on the uniform bound on ||wj||l°°(s(_r)) such that for any 

(6.3) \\diWejUj\\LP{B{B.) + \\d2^ejU3\\LP{B{R) 

< W^lj"^ e,Uj\\LP{B{R) + 11(1 + \Uj\)d2Ve,Uj\\LP{BiR) 

^ \\'^j\\w^'''{B{R)) +C'l||Wj|lw/^P(i3(i?)) +C2<C, 

for a new constant C > independent of j. The weak regularity of u and the weak Sobolev convergence 
follow in a standard fashion. 

Next we address the PDE: Since for every j the approximating solution Uj is of class C°° then we can use 
the non divergence form of the equation 

2 

2J 0,hk{^ ]Uj)Xh^]Xk,jUj = 0. 
h.k=l 

Here 

ah,kiyjUj) -* a/i,fc(Vow) = ShiSki in L^, 
while 

(6.4) ^ijWj -^ Xu, X2,jUj ^ 

as j — > +00 weakly in Wi^^ (O) . Hence letting j go to 00 in the non divergence form equation we conclude 

X^u = 
in the Sobolev sense. D 



An analogous result holds for higher order derivatives: 

Proposition 6.2. For every k ^ N for every p > 1 and for every multiindex I of length k, the sequence 
iyl.Uj) is bounded in W^£;'^oc(^)- Moreover 

X'^jUj -^ X'^u, and X^jUj -> 

weakly in W-^'^{fl) as j —> 00. We will express this convergence in the notation 

^i^.Uj — > DqU as j -^ +00, weakly in W^'^{fl). 
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Proof. Arguing as in the previous result and using (I5.14p . and ()6.ip we deduce that for every ball B{R) CC fi 
there exist Ci, C2, C independent of j, such that 

\\d2Dl^U,\\L.iBiR)) <Cl J2 \\De,Vj\\L.iBiR)) + C^ < C , 

\J\<\I\ 

\\9xDI^Uj\\lp{B(R)) < \\Xl,jDlMj\\LPiB{R)) + \\ujd2Dl^Uj\\LP{B(R)) < C, 

SO that (V^.Uj) is bounded in W^'^^^{^) for every p > 1, and every niulti- index /. Let us prove that the 
weak limit is DqU. If |/| = 1 the assertion is true by (|6.2p . If / is a multi-index such that |/| = k, we can 
assume by simplicity that / = (1, /'), where |/'| = k — 1. We can also assume by inductive hypothesis that 

d2Uj ~* d2U as J — > 00 weakly in L^ifl) 

Uj — > u as j ~> 00 in Lf^^(fl) 

y g-Uj -^ Dq u as j ^ 00 in Lf^^{Q). 
Then integrating by parts 

lim / W^.Ujcf) ^ — lim / W^.UjXjCJ) ^ / d2UjW^.Uj(j> = 

]-*oo J J 3-*ooJ ' J J 

Dq uX(j) — / d2uDQ u(j), 
and this ensures the weak convergence of {DI.Uj) to D^u. D 

Remark 6.3. In view of the Ascoli-Arzeld theorem and the Morrey-Sobolev embedding one has convergence 
Xi aUj -^ X^u in the C" norm on compact subsets of Q. for all a G (0, 1). 

We can now prove the main regularity properties of the limit function u: 

Proposition 6.4. For every k, and for every p > 1 the function z = X^u belongs to W^^i^J^Vt) and it is an 
a.e. solution of 

(6.5) X^z^O in n. 
In particular 

(6.6) X'^M e cr^^m 

for every < a < 1. 

Proof. Since u is a vanishing viscosity solution of X^u = in O, then Proposition 16.21 implies X^u G 
^B loc(^) ^*-"^ all p > 1. As X^u = a.e. in fl, then a simple iteration shows that all the derivatives 
X'^X'^u vanish a.e. in fl. The Holder regularity (|6.6|) follows from the classical Morrey-Sobolev embedding 
theorem. D 

We can now give a new pointwise definition of derivative in the direction of vector fields Xi and A2. 

Definition 6.5. Let X be a Lipschitz vector field on Q and let S,o (z Q and j be a solution to problem 
7' = X(7), 7(0)=^. 

We say that a function f e C;"^(r2), with a e]0, 1[, has Lie-derivative in the direction of the vector field 
X in ^0 if there exists 

^(/O7)|.=o, 
and we will denote its value by Xf{£^o). 
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If the weak derivative of a function / is sufHciently regular, then the two notions of derivatives coincide. 
For the proof of the following result see [131 Remark 5.6]. 

Proposition 6.6. // / G C"^^{n) for some a e]0, 1[ and its weak derivatives Xf e C^^^{Q),d2f G L^^^{Q) 
with p > 1/a, then for a// ^ G il the Lie-derivatives X f(^) exist and coincide with the weak ones. 

We are now ready to prove the result concerning the foliation 
Proof of Corollary 11.61 The equation 7' = X{j) has an unique solution, of the form 

7(2;) = {x,y{x)), 

where y'{x) = u{x,y{x)). In view of the regularity of u and of the previous proposition then y"{x) = 
Xu{x, y{x)), and y"'{x) = X'^u{x, y{x)) = 0. This shows that 7 is a polynomial of order 2 and concludes the 
proof. n 
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